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RANS Modeling of Benchmark Shockwave / Boundary 

Layer Interaction Experiments 


Abstract: This presentation summarizes the computations of a set of shock wave / 
turbulent boundary layer interaction (SWTBLI) test cases using the Wind-US code, as 
part of the 2010 American Institute of Aeronautics and Astronautics (AIAA) shock / 
boundary layer interaction workshop. The experiments involve supersonic flows in 
wind tunnels with a shock generator that directs an oblique shock wave toward the 
boundary layer along one of the walls of the wind tunnel. The Wind-US calculations 
utilized structured grid computations performed in Reynolds-averaged Navier-Stokes 
mode. Three turbulence models were investigated: the Spalart-Allmaras one- 
equation model, the Menter Shear Stress Transport k-co two-equation model, and an 
explicit algebraic stress k-co formulation. Effects of grid resolution and upwinding 
scheme were also considered. The results from the CFD calculations are compared 
to particle image velocimetry (PIV) data from the experiments. As expected, 
turbulence model effects dominated the accuracy of the solutions with upwinding 
scheme selection indicating minimal effects. 


1 


National Aeronautics and Space Administration - Glenn Research Center 



RANS Modeling of Benchmark 
Shockwave / Boundary Layer 
Interaction Experiments 

Nick Georgiadis, Manan Vyas, and Dennis Yoder 

RTE/Inlet and Nozzle Branch 
NASA Glenn Research Center 
Cleveland, OH, USA 


2 


National Aeronautics and Space Administration - Glenn Research Center 


Outline 

• U. Michigan Case 1 - 7.75° Shock Generator, Mach 2.75 

- Grid and solution details with SST turbulence model 

- Comparisons of CFD solutions with experiment using SST, SA, 
ASM turbulence models 

- Upwind scheme comparisons with SST 

- Grid resolution comparisons with SST 

• U. Michigan Case 1-3 (7.75°, 10°, and 12°) U-velocity 
variations, all Mach 2.75 

• UFAST*, 8° Shock Generator, Mach 2.25 

- Grid and solution details with SST turbulence model 

- Comparisons of CFD solutions with experiment using SST, SA, 
ASM turbulence models 


*Unsteady effects in shock wove induced separation, Institut Universitaire des Systemes 
Thermiques Industriels 3 
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Flow Solver Specifics 

• Wind-US, version 2.0 

• Structured grids utilized, grids built with Gridgen 

• Turbulence Models: 

- SST* (Menter Shear Stress Transport) 

- SA (Spalart Allmaras) 

- ASM (built upon Menter BSL k-co) 

• Upwind Schemes (all 2 nd order): Roe*, Van-Leer, HLLC 


* items shown in bold are default when not otherwise noted. 
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Grid (1 of 3) 





• 8.4 M grid points, 16 structured zones 

• point-to-point connectivity in shock generator region 

• y+ = 1 .2 to 1 .5 for no-slip walls (calculated using 
freestream conditions and “average” C f =0.0025) 

•Ax+ ~ 100 at key axial stations (i.e. changes in 
geometry) 

• Grids smooth (no abrupt changes in packing) in shock 
generator region, within & across zones 
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UM Cases 1-3 Grid (2 of 3) 
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UM Cases 1-3 Grid (3 of 3) 



Entrance plane to shock 
generator region: 



National Aeronautics and Space Administration - Glenn Research Center 



U. Michigan Case 1 - 7.75° Shock Generator, Mach 2.75 

Comparisons of CFD solutions with experiment 
using SST, SA, ASM turbulence models 
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U Velocities Along Symmetry Plane 

(7.75 degree wedge) 
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Eddy Viscosity Along Symmetry Plane 

(7.75 degree wedge) 
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* Contours shown are ix t /\i ref and are clipped on high end. 
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Eddy Viscosity Along Symmetry Plane 
(7.75 degree wedge) Old 
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* Contours shown are fi t /fi ref and are clipped on high end. 
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Eddy Viscosity Along Symmetry Plane 

(10 degree wedge) NEW 
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* Contours shown are fi t /fi ref and are clipped on high end. 
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Eddy Viscosity Along Symmetry Plane 

(12 degree wedge) Old 



Y 
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* Contours shown are fi t /fi ref and are clipped on high end. 
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Eddy Viscosity Along Symmetry Plane 

(10 degree wedge) Old 



Y 
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* Contours shown are fi t /fi ref and are clipped on high end. 
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U Velocity Contours 



Experiment: 




SA: 


ASM: 
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V Velocity Contours 
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dV/dX Contours 
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U Velocity Profiles 

x = 20.76 mm x = 30.76 mm 
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U Velocity Profiles 


x = 38.76 mm 


x = 53.76 mm 





19 


y (mm) 


National Aeronautics and Space Administration - Glenn Research Center 

V Velocity Profiles 

x = 20.76 mm x = 30.76 mm 
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V Velocity Profiles 

x = 38.76 mm x = 53.76 mm 
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U. Michigan Case 1 - 7.75° Shock Generator, Mach 2.75 

Upwind scheme comparisons with SST 
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U Velocity Profiles 



(Effect of Upwind Scheme - all SST Turbulence Model) 


x = 20.76 mm x = 53.76 mm 
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V Velocity Profiles 

(Effect of Upwind Scheme - all SST Turbulence Model) 

x = 20.76 mm x = 53.76 mm 
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U. Michigan Case 1 - 7.75° Shock Generator, Mach 2.75 

Grid resolution comparisons with SST, Roe Scheme 

• Fine grid = 8.4 M points 

• Medium grid = 1.1 M points (every other point utilized in each 

direction). 

• Coarse (every 4 th point in each direction) also available, but not 
shown in comparisons here. 
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U Velocity Profiles 

x = 20.76 mm x = 53.76 mm 
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V Velocity Profiles 


x = 20.76 mm 


x = 53.76 mm 
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U. Michigan Case 1-3 (7.75°, 10°, and 12°) Mach 2.75 

U-velocity variations 
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U Velocities Along Symmetry Plane 

(7.75 degree wedge) 
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U Velocities Along Symmetry Plane 

(10 degree wedge) 
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U Velocities Along Symmetry Plane 
(10 degree wedge Manan) 
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U Velocities Along Symmetry Plane 

(12 degree wedge) 
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UFAST 8° Shock Generator, Mach 2.25 

- Grid and solution details with SST turbulence model 

- Comparisons of CFD solutions with experiment using SST, SA, 
ASM turbulence models 
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UFAST Grid (1 of 2) 





• 7.3 M grid points, 10 structured zones 

• point-to-point connectivity in shock generator region 

• y+ = 1 .0 for no-slip walls (calculated using freestream 
conditions and “average” C f =0.0025) 

•Ax+ ~ 60 at key axial stations (i.e. changes in 
geometry) 

• Grids smooth (no abrupt changes in packing) in shock 
generator region, within & across zones 


National Aeronautics and Space Administration - Glenn Research Center 

UFAST Grid (2 of 2) 



Zoom on shock 
generator region 
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U Velocities Along Symmetry Plane 
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Eddy Viscosity Along Symmetry Plane 
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* Contours shown are p t /p ref and are clipped on high end. 
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U Velocity Contours 


Experiment SST: 




SA: 


ASM: 
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U Velocity Profiles 

x = 280 mm x = 300 mm 
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U Velocity Profiles 


x = 320 mm 
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Summary of Cases Completed with Wind-US 


• Solutions obtained with SST and SA for all 4 cases 
(3 U. Michigan cases & UFAST Case) on medium and 
fine grid levels. 

• ASM solutions obtained for Michigan case 1 (7.75°) & 
UFAST case, medium and fine grid levels. 

• Numerical scheme investigations for U. Michigan case 1 . 
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Normalized Wall Pressure, P/P(ref) 
Experimental P(ref) @ Isolator Entrance 
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UVA- Direct Connect Sc ram jet 



STATIC PRESSURE PROFILE 
Clean Air, Fuel-off, Medium Grid 
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Conclusions 



• Turbulence model effects are dominant, as expected. 

• SST performs “best” however not great agreement. 
Further, our other experiences show SST failing for 
cases with more severe shocks & separations (i.e. 
isolators, overexpanded nozzles). 

• Minimal effects of numerical scheme. 

• Grid density investigations examined sequenced grids. 

• Grid density investigations only examined sequenced 
grid. A more comprehensive study would target specific 
grid clustering to determine shock resolution needs. 
Current fine grid seems adequate, but finer grid needed 
to identify grid convergence. 
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